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As far as the laws of mathematics refer to 
reality, they are not certain; and as far as 
they are certain, they do not refer to reality. 

A. Einstein 



Abstract 

An interesting example of the deep interrelation between Physics and Mathematics is 
obtained when trying to impose mathematical boundary conditions on physical quantum 
fields. This procedure has recently been re-examined with care. Comments on that and 
previous analysis are here provided, together with considerations on the results of the purely 
mathematical zeta-function method, in an attempt at clarifying the issue. Hadamard reg- 
ularization is invoked in order to fill the gap between the infinities appearing in the QFT 
renormalized results and the finite values obtained in the literature with other procedures. 
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1 Introduction 



The question, most beautifully phrased by Eugene Wigner as that of the unreasonable effec- 
tiveness of mathematics in the natural sciences [1] is an old and intriguing one. It certainly 
goes back to Pythagoras and his school ( "all things are numbers"), even probably to the 
sumerians, and maybe to more ancient cultures, which left no trace. In more recent times, 
I. Kant and A. Einstein, among others, contributed to this idea with extremely beauti- 
ful and profound reflections. Also, mathematical simplicity, and beauty, have remained for 
many years, and to many important theoretical physicists, crucial ingredients when having 
to choose among different plausible possibilities. 

A crude example of unreasonable effectiveness is provided by the regularization pro- 
cedures in quantum field theory based upon analytic continuation in the complex plane 
(dimensional, heat-kernel, zeta-function regularization, and the like). That one obtains a 
physical, experimentally measurable, and extremely precise result after these weird mathe- 
matical manipulations is, if not unreasonable, certainly very mysterious. No wonder that for 
more one highly honorable physicist, the taming of the infinities that plague the calculations 
of modern quantum physics and the ensuing regularization/renormalization issue, were al- 
ways unacceptable, illegal practices. Those methods are presently full justified and blessed 
with Nobel Prizes, but more because of the many and very precise experimental checkouts 
(the effectiveness) than for their intrinsic reasonableness. 

As this note is meant not just for specialists, a simple example may be clarifying. Consider 
the calculation of the zero point energy (vacuum to vacuum transition, also called Casimir 
energy [2]) corresponding to a quantum operator, H, with eigenvalues A n , 

E =<0\H\0>=W X m (1) 

1 n 

where the sum over n is a sum over the whole spectrum (which will involve, in general, several 
continuum and several discrete indices.) The last appear typically when compactifying 
the space coordinates (much in the same way as time compactification gives rise to finite- 
temperature field theory). In turns out, in fact, that only in very special cases will this sum 
yield a convergent series. Generically one has to deal with a divergent sum, which is then 
regularized by different appropriate means. The zeta-function method [3] — which stands on 
solid and flourishing mathematical grounds [4] — will interprete this sum as being the value 
obtained for the zeta function of the operator H, namely 

U*) = EV, (2) 

n 

at the point s = — 1 (we set ft = c = 1). It turns out generically that (h(s) is well defined 
as an absolutely convergent series for Re s > a (a is a certain abscissa of convergence), 
and that it can be analytically continued (in a perfectly defined way) to the whole complex 
plane, with the possible appearance of poles as only singularities. If, as it often happens, 
Cft(s) has no pole at s = — 1 we are done; if it hits a pole there, then further elaboration 
of the method is necessary. That the mathematical result one gets through this process 
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of analytical continuation in the complex plane, does coincide with the experimental result, 
constitutes the specific example of unreasonable effectiveness of mathematics we are referring 
to. 

As one should expect, things do not turn out to be so simple and clear-cut in practice. 
Actually, one cannot assign a meaning to the absolute value of the zero-point energy, and 
any physically measurable effect is to be worked out as an energy difference between two 
situations, such as a quantum field in curved space as compared with the same field in flat 
space, or a field satisfying boundary conditions on some surface as compared with the same 
field in the absence of the surface, etc. This energy difference is the Casimir energy: 



And it is here where the problem appears. Imposing mathematical boundary conditions 
on physical quantum fields turns out to be a highly non-trivial act. This was discussed in 
much detail in a paper by Deutsch and Candelas a quarter of a century ago [5] . These authors 
quantized electromagnetic and scalar fields in the region near an arbitrary smooth boundary, 
and calculated the renormalized vacuum expectation value of the stress-energy tensor, to find 
that the energy density diverges as the boundary is approached. Therefore, regularization 
and renormalization did not seem to cure the problem with infinities in this case and an 
infinite physical energy was obtained if the mathematical boundary conditions were to be 
fulfilled. However, the authors argued that, in nature, surfaces have non-zero depth, and 
its value could be taken as a handy (dimensional) cutoff in order to regularize the infinities. 
This approach will be recovered later in the present paper. Just two years after Deutsch and 
Candelas' work, Kurt Symanzik carried out a rigorous analysis of quantum field theory in the 
presence of boundaries [6] . Prescribing the value of the quantum field on a boundary means 
using the Schrodinger representation, and Symanzik was able to show rigorously that such 
representation exists to all orders in the perturbative expansion. He showed also that the field 
operator being diagonalized in a smooth hypersurface differs from the usual renormalized 
one by a factor that diverges logarithmically when the distance to the hypersurface goes to 
zero. This requires a precise limiting procedure, and also point splitting, to be applied. In 
any case, the issue was proven to be perfectly meaningful within the domains of renormalized 
quantum field theory. One should note that in this case the boundary conditions and the 
hypersurfaces themselves were always treated at the pure mathematical level (zero depth) 
by using delta functions. 

More recently, a new approach to the problem has been postulated [7] which uses elements 
of the two methods above. Boundary conditions on a field, 0, are enforced on a surface, S, 
by introducing a scalar potential, a, of Gaussian shape living on the surface (and near it). In 
the limit when the Gaussian becomes a delta function, the boundary conditions (Dirichlet, 
in this case) are enforced, in that the delta-shaped potential kills all the modes of at the 
surface. For the rest, the quantum system undergoes a full-fledged quantum field theory 
renormalization, as in the case of Symanzik's approach. The results obtained confirm those 
of [5], in the several models studied albeit they do not seem to agree with those of [6]. 




(3) 
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Such results are also in clear contradiction with the ones quoted in the usual textbooks 
and review articles dealing with the Casimir effect [8], where no infinite energy density when 
approaching any of the Casimir plates has been reported. 

2 A zeta-function approach 

Too often has it been argued that sophisticated regularization methods, as for instance 
the zeta-function regularization procedure, get rid of the infinities in a rather obscure way 
(through analytic continuation, as previously described), so that, contrary to what happens 
when a cut off is introduced, one cannot keep trace of the infinites, which are cleared up 
without control, leading sometimes to erroneous results. 

There is no way to refute an statement of this kind rigorously, but it should be pointed 
out that in many ocasions (if not always) the discrepancies between the result obtained 
by using the zeta procedure and other, say cut-off like, approaches have been proven to 
emerge from a misuse of zeta regularization, and not to stem from the method itself. When 
employed properly, the correct results have been recovered (for a good number of examples, 
see [3, 4, 9, 10]). 

Take the most simple example of a scalar field in one dimension, (f>(x), with a boundary 
condition (BC) of Dirichlet type imposed at a point (say x — 0), e.g. 0(0) = 0. We would like 
to calculate the Casimir energy corresponding to this configuration, that is, the difference 
between the zero point energy corresponding to this field when the boundary condition is 
enforced, and the zero point energy in the absence of any boundary condition at all. Taken 
at face value, both energies are infinite. The regularized difference may still be infinite when 
the BC point is approached (this is the result in [7]) or might turn out to be finite (even 
zero, which is the result given in some standard books on the subject). 

Let us try to understand this discrepancy. According to Eq. (1), we have to add up all 
energy modes. For the mode with energy u, the field equation reduces to: 

-4>"(x) + m 2 (j)(x) = oo 2 (j)(x). (4) 

In the absence of a BC, the solutions to the field equation can be labeled by k = +\Juj 2 — m 2 > 
0, as 

<P k (x) = Ae lkx + Be~ ikx , (5) 

with A, B arbitrary complex (for the general complex), or as 

4>k{x) = asin(A;:r) + bcos(kx), (6) 

with a, b arbitrary real (for the general real solution). Now, when the mathematical BC of 
Dirichlet type, 0(0) = 0, is imposed, this does not influence at all the eigenvalues, k, which 
remain exactly the same (as stressed in the literature). However, the number of solutions 
corresponding to each eigenvalue is reduced by one half to: 

0l D) (x) = A(e fc -e-^), (7) 
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with A arbitrary complex (complex solution), and 

4 >( k > \ x ) — asin(kx), (8) 

with a arbitray real (real solution). In other words, the energy spectrum (for omega) that 
we obtain in both cases is the same, a continuous spectrum 



Vm 2 + k 2 , (9) 



but the number of eigenstates corresponding to a given eigenvalue is twice as big in the 
absence of the BC. 2 

Of course these considerations are elementary, but they seem to have been put aside 
sometimes. They are crucial when trying to calculate (or just to give sense to) the Casimir 
energy density and force. More to the point, just in the same way as the traces of the two 
matrices Mi = diag (a, (3) and M2 = diag (a, a, /3, (3) are not equal in spite of having "the 
same spectrum a, (3," it also turns out that, in the problem under discussion, the traces of 
the Hamiltonian with and without the Diriclet BC imposed yield different results, both of 
them divergent, namely 

roo 

tr H = 2 tr H BC = 2 dk Vm 2 + k 2 . (10) 

Jo 

By using the zeta function, we define 

( BC (s):= rdK^ + J)-', »:=-, (11) 
Jo [l 

with n a regularization parameter with dimensions of mass. 3 We get 

CBc( S )= v¥r 2 ( ;- 1/2) (^)^-, (12) 

and consequently, 

t?H BC = 1 -Cbc{s = -1/2) 



2 

m 



4^ 



I m 2 

-jiz + 1 - 7 - log — - *(-l/2) + 0(s + 1/2) 



■ (13) 

-1/2 



As is obvious, this divergence is not cured when taking the difference of the two traces, Eq. 
(3), in order to obtain the Casimir energy: 



E c /v = E* c /^ - E /» = -E^'/^ = V -^^. (14) 



2 To understand this point even better (by making recourse to what is learned in the maths classes at 
high school), consider the fact that further, by imposing Cauchy BC: (f)(0) — 0,(f>'(0) — 0, the eigenvalues 
remain the same, but for any k the family of eigenfunctions shrinks to just the trivial one: <j>k{x) = 0, VTc 
(the Cauchy problem is an initial value problem, which completely determines the solution). 

3 Always necessary in zeta regularization, since the complex powers of the spectrum of a (pseudo-) dif- 
ferential operator can only be defined, physically, if the operator is rendered dimensionless, what is done by 
introducing this parameter. That is also an important issue, which is sometimes overlooked. 
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Wejust hit the pole of the zeta function, in this case. 

How is this infinite to be interpreted? What is its origin? Just by making recourse to the 
pure mathematical theory (reine Mathematik) , we already get a perfect description of what 
happens and understand well where does this infinite energy 4 come from. It clearly originates 
from the fact that imposing the boundary condition has reduced to one-half the family of 
eigenf unctions corresponding to any of the eigenvalues which constitute the spectrum of 
the operator. And we can also advance that, since this dramatic reduction of the family 
of eigenfunctions takes place precisely at the point where the BC is imposed, the physical 
divergence (infinite energy) will originate right there, and nowhere else. 

While the analysis above cannot be taken as a substitute for the actual modelization of 
Jaffe et al. [7] — where the BC is explicitly enforced through the introduction of an auxiliary, 
localized field, which probes what happens at the boundary in a much more precise way — 
it certainly shows that pure mathematical considerations, which include the use of analytic 
continuation by means of the zeta function, are in no way blind to the infinites of the physical 
model and do not produce misleading results, when the mathematics are used properly. And 
it is very remarkable to realize how close the mathematical description of the appearance of 
an infinite contribution is to the one provided by the physical realization [7]. 



3 The case of two-point Dirichlet boundary conditions 

A similar analysis can be done for the case of a two-point Dirichlet BC: 

0(a) = 0, 0(-a) = 0. (15) 

Straightforward algebra shows, in this situation, that the eigenvalues k are quantized, as 
k — 7r/ (2a), so that: 



^ = \jm 2 + —, ^ = 0,1,2,... (16) 

The family of eigenfunctions corresponding to a given eigenvalue, u>e, is of continuous dimen- 
sion 1, exactly as in the former case of a one-point Dirichlet BC, namely, 

<f>t(x) = fosin - a^j , (17) 

where b is an arbitrary, real parameter. 5 To repeat, the act of imposing Dirichlet BC on 
two points has the effect of discretizing the spectrum (as is well known) but there is no 
further shrinking in the number of eigenfunctions corresponding to a given (now discrete) 
eigenvalue. 

The calculation of the Casimir energy, by means of the zeta function, proceeds in this case 
as follows [3, 4, 9, 10]. To begin with, it may be interesting to recall that the zeta-'measure' 

4 In mathematical terms, this infinite value for the trace of the Hamiltonian operator. 
5 The contribution of the zero-mode (I = 0) is controverted, but we are not going to discuss this issue 
here (see e.g. [11] and references therein.) 
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of the continuum equals twice the zeta-'measure' of the discrete. In fact, just consider the 
following regularizations: 



n=l n=l 

and 

,1-s 



= Kr(0) = -$, (18) 



poo poo 

/ dk= dk(k + fx)- 
Jfi Jo 



: £ 

s=0 S - 1 



= -ii. (19) 

s=0 



The result is, as advanced, that the ( 'measure' of the discrete is half that of the continuum. 

The trace of the Hamiltonian corresponding to the quantum system with the BC imposed, 
in the massive case, is obtained by means of the zeta function 



t BC (s) ■= E 7T + 



00 'm 2 n 2 f \ s 



f^[\n 2 4:fi 2 a 2 J 
(»y\ 1 T(s - 1/2) am 2n° /2am\V 2 +* - 1/2 

U -2 + —[x^^ + my—) S s - 1/2( 



anm 



, (20) 



being K v a modified Bessel function of the third kind (or MacDonald's function). Thus, for 
the zero point energy of the system with two-point Dirichlet BC, we get 

tr HBO/? = -(bc(s = -1/2) = --^jf- " ^ E -Ki(2*nm/ri, (21) 

where \i is, in this case, /i := n/(2a) (a fixes the mass scale in a natural way here). As in 
the previous example, we finally obtain an infinite value for the Casimir energy, namely 

FY — 1 )m 2 777 00 1 

Ed ii = E* c lv - E /v = ~ £^ E -Kifrmm/l*)- (22) 

It is, therefore, not true that regularization methods using analytical continuation (and, 
in particular, the zeta approach) are unable to see the infinite energy that is generated on 
the boundary-conditions surface [5, 6, 7] (see Eq. (34) later). The reason is still the same as 
in the previous example: imposing a two-point Dirichlet BC amounts again to halving the 
family of eigenfunctions which correspond to any given eigenvalue (all are discrete, in the 
present case, but this makes no difference). In physical terms, that means having to apply 
an infinite amount of energy on the BC sites, in order to enforce the BC. In absolute analogy, 
from the mathematical viewpoint, halving the family of eigenfunctions immediately results 
in the appearance of an infinite contribution, under the form of a pole of the zeta function. 

The reason why these infinities (the one here and that in the previous section) do not 
usually show up in the literature on the Casimir effect is probably because the textbooks on 
the subject focus towards the calculation of the Casimir force, which is obtained by taking 
minus the derivative of the energy with respect to the plate (or point) separation (here 
w.r.t. 2a). Since the infinite terms do not depend on a, they do not contribute to the force (as 
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is recognized explicitly in [7]). However, some erroneous statements have indeed appeared 
in the above mentioned classical references, stemming from the lack of recognition of the 
catastrophical implications of the act of halving the number of eigenf unctions, when imposing 
the BC. The persistence of the eigenvalues of the spectrum was probably misleading. We 
hope to have clarified this issue here. 

4 How to deal with the infinities 

Here, the infinite contributions have shown up at the regularization level, but a more careful 
study [7] was able to prove that they do not disappear even after renormalizing in a proper 
way. The important question is now: are these infinities physical! Will they be observed as 
a manifestation of a very large energy pressure when approaching the BC surface in a lab 
experiment? No doubt such questions will be best answered in that way, e.g. experimentally. 
If, on the contrary, that sort of large pressures fails to manifest itself, this might be a clear 
indication of the need for an additional regularization prescription. In principle, this seems 
to be forbidden by standard renormalization theory, since the procedure has been already 
carried out to the very end: there remains no additional physical quantity which could 
possibly absorb the divergences (see [7]). 

In any case, there are circumstances — both in physics and in mathematics — where cer- 
tain 'non-orthodox' regularization methods have been employed with promising success. In 
particular, Hadamard regularization in higher-post-Newtonian general relativity [12] and also 
in recent variants of axiomatic and constructive quantum field theory [13]. Among mathe- 
maticians, Hadamard regularization is nowadays a rather standard technique in order to deal 
with singular differential and integral equations with boundary conditions, both analytically 
and numerically (for a sample of references see [14]). Indeed, Hadamard regularization is a 
well-established procedure in order to give sense to infinite integrals. It is not to be found in 
the classical books on infinite calculus by Hardy or Knopp; it was L. Schwartz [15] who pop- 
ularized it, rescuing Hadamard's original papers. Nowadays, Hadamard convergence is one 
of the cornerstones in the rigorous formulation of QFT through micro-localization, which on 
its turn is considered by specialists to be the most important step towards the understanding 
of linear PDEs since the invention of distributions (for a beautiful, updated treatment of 
Hadamard's regularization see [16]). 

Let us briefly recall this formulation. Consider a function, g(x), expandable as 

g(x) = j:-^- r + h(x), (23) 

with Xj complex in general and h(x) a regular function. Then, it is immediate that 

[" dxg(x) = P(l/e) + H(e), (24) 

J a+e 

being P a polynomial and H(0) finite. If the Xj N, then one defines the Hadamard 
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regularized integral as 



r° r° a- 

=h dxg(x) := / h(x) dx — ^ T — J ~~r(^ ^ a ) 

Ja Ja „■_-. \j — 1 



1-A,- 



Alternatively, one may define, for a N, p < a < p + 1, and e C[-i,i] 



(25) 



J ' r(-a)T-i (i-tY+^ 



V(-a)J-i (1-t) 

to obtain, after some steps, 

/ (J) (-1) 

^r0' + i-a)2«-j ' r(p 



(26) 



(27) 



where the last integral is at worst improper (Cauchy's principal part). If Ai = 1, then the 
result is a\ In (6 — a), instead. If Ai = p e N, calling 



H p (f;x) := £ 



dt 



/(*) 



(t-x) 



p+i' 



\x\ < 1, 



(28) 



we get 
H p (f;x) 



p fU)( T ) 



dt 



(t - x)p +1 



dt 



3=0 J - 

where the first term is regular and the second one can be easily reduced to 

1 d p ' j r 1 dt 



(t - x)p +1 -i ' 



(29) 



{p — j)\ dx p ~i 



f l dt 

> i-i t-x' 



(30) 



being the last integral, as before, a Cauchy PP. 

An alternative form of Hadamard's regularization, which is more fashionable for phys- 
ical applications (as is apparent from the expression itself) is the following [12]. For the 
case of two singularities, at xi, x 2 , after excising from space two little balls around them, 
R 3 \ (B ri (xi) U B r2 (x 2 )), with B ri (xi) fl B r2 (x 2 ) = 0, one defines the regularized integral as 
being the limit 



j d 3 x F(x) := FP Q ,^ / d 3 x F(x), 



(31) 



where S\ and s 2 are two (dimensionfull) regularization parameters [12]. This is the version 
that will be employed in what follows. 



5 Hadamard regularization of the Casimir effect 

We now use Hadamard's regularization as an additional tool in order to make sense of the 
infinite expressions encountered in the boundary value problems considered before. As it 
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turns out from a detailed analysis of the results in [7] (which we shall not repeat here, 
for conciseness), the basic integrals which produce infinities, in the one-dimensional and 
two-dimensional cases there considered, are the following. 

In one dimension, with Dirichlet BC imposed at one (x = 0) and two (x = ±a) points, 
respectively, by means of a delta-background of strength A (see [7]), one encounters the two 
divergent integrals: 



E 1 (X,m) = 
E 2 (a,X,m) = 



1 

2tt 



dt 



1 r c 

2lT Jm 



dt 




A 



(32) 
(33) 



Using Hadamard's regularization, as described before, we obtain for the first one, Eq. (32), 

f • < 34 > 



Ei(m)= ^-(l-ln- 
471 \ m 



+ 



A^OO 



where the first term is the singular part when the limit A 
—which is Hadamard's finite part — yields in this case 



m 



oo is taken, and the second 



(35) 



Such result is coinciding with the classical one (0, for m — 0). Note in particular, that the 
further mm divergence as m — > oo is hidden in the A— divergent part, and that behavior 
does explain why the classical results which are obtained using hard Dirichlet BC — what 
corresponds as we just prove here to the Hadamard's regularized part — cannot see anything 
of it. 

In the case of a two-point boundary at x — ±a (separation 2a), Eq. (33), we get a similar 
Eq. (34) but now the regularized integral is as follows. For the massless case, we obtain 



i 



48a' 



(36) 



which is the regularized result to be found in the classical books. In the massive case, m^fl, 
after some additional work the following fast convergent series turns up [cf. Eq. (22)] 



i 



m °° 1 
-^Y.-^dakm). 

fc=l 



(37) 



Thus Eq. (34) yields strictly the same result (22) that was already obtained by imposing 
the Dirichlet BC ab initio. What has now been gained is a more clear identification of the 
singular part, in terms of the strength of the delta potential at the boundary. This will be 
the general conclusion, common to all the other cases here considered. 
Correspondingly, for the Casimir force we obtain the finite values 6 

7T 



^2 (a) 



96a 2 



(38) 



6 Note that the force F(a) is here given by minus the derivative of the total energy E{a) w.r.t. 2a, since 
this is the distance between the two Dirichlet points {not a.) 
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in the massless case, and in the massive one 



F 2 (a,m) = 2^ 



2 oo 



7T 



k=l 



K {Aakm) + 



Aakm 



K x [Aakm) 



(39) 



Those expressions coincide with the ones derived in the above mentioned textbooks on the 
Casimir effect, and reproduced before by using the zeta function method (just take minus 
the derivative of Eq. (34) w.r.t. 2a). 

The two dimensional case turns out to be more singular [7] — in part just for dimensional 
reasons — and requires additional wishful thinking in order to deal with the circular delta 
function sitting on the circumference where the Dirichlet BC are imposed. Here one encoun- 
ters the basic singular integral, for the term contributing to the second Born approximation 
(we use the same notation as in [7]), 



[•OO 

&(p) = / dr rJ (pr)a(r), air) =6Aexp 
Jo 

with J a Bessel function of the first kind, and 

roo uj^o 

/ drair) — A, air) X5(r — a). 

Jo 



(r - a)' 
2uj 2 



(40) 



(41) 



Hadamard's regularization yields, in this case (the r's replacing the cr's in the regularized 
version) , 



r(r,p) = c\(rp + 1) ^exp 



(r - a)' 
2u 2 



LU^O 



A 5(r — a), 



(42) 



with p a (dimensionfull) regularization parameter, being the constant c given by 

-i f°° , , . (r - a) 2 

c = arr exp — 

Jo ' 2ur 

which exists and is perfectly finite; in particular, c~ l {u — .l,a — 1) — .25. Then, 

roo 

fip) = 2n drr Jo(pr) r(r,p) = 2n\a(ap + l)~ w ' 2 J (ap) 
Jo 

It turns out that for the Casimir energy we get in this case (notation as in [7]) 
X 2 a 2 r°° 



(43) 



(44) 



8 Jo 
\ 2 a 2 < 



dp(ap+l) w Jo(ap) arctan(p/2m) 



LO^O 



8 2uj 



1 7+3m2 
h - — h 4m 



2a 



7 -= [1 — ln(am)] h(Aa 2 m 2 ) 

V 71 " 



(45) 



where h(z) := 2-^3 ((1/2, 1/2); (1, 1, 3/2); z) and 7 is the Euler-Mascheroni constant; in par- 
ticular, for instance h(l) = 1.186711, what is quite a nice value. Recall also that u is the 
width of the Gaussian 5, which is the very physical parameter considered in [5]. When this 
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width tends to zero an infinite energy appears (the width controls the formation of the pole). 
The rest of the result is the Hadamard regularization of the integral, e.g. 7 



Again, the finite part reverts to the results obtained in the literature with Dirichlet BC ab 

initio. 

To summarize, it has been here proven — in some particular but rather non-trivial and 
representative examples — that the finite results derived through the use of Hadamard's 
regularization exactly coincide with the values obtained using the more classical, less full- 
fledged methods to be found in the literature on the Casimir effect. Moreover, Hadamard's 
prescription is able to separate and identify the singularities as physically meaningful cut- 
offs. Although the validity of this additional regularization is at present questionable, the 
fact that it bridges the two approaches is already remarkable, maybe again a manifestation 
of the unreasonable effectiveness of mathematics. 
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